How fast could a proto-pulsar rotate? 
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O ■ ABSTRACT 

c : 

| ^ . According to two estimated relations between the initial period and the 

. dynamo-generated magnetic dipole field of pulsars, we calculate the statistical 

distributions of pulsar initial periods. It is found that proto-pulsars are very 
i . likely to have rotation periods between 20 and 30 ms, and that most of the 

^ , pulsars rotate initially at a period < 60 ms. 
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\q ■ Subject headings: pulsars — neutron stars — magnetic fields 

o : 
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^3 ■ 1. Introduction 

Pulsars provide us an unusual physical condition to increase our knowledge of the nature. 
As yet, one of the essential parameters, the initial periods Pq of proto-pulsars, 1 which may 
reveal precious information of dynamical supernova process, are poorly known. Actually 
there are 3 efforts to estimate Pq in the literatures. 1, The period Po can be found via 
T = P/[(n - 1)P][1 - (Po/P)™- 1 ] if the age T and the braking index n = QQ/Q 2 (Q = 2vr/P 
the angular velocity of rotation) are measured (e.g., Kaspi et al. 1994), where P is the 
rotation period observed. This method needs constant braking indices, which are likely to 
vary with time (Xu & Qiao 2001). 2, Monte Carlo simulation is used to study the pulsar 
"current" in the magnetic field - period diagram (e.g., Lorimer et al. 1993); the initial period 
of "injected" pulsars into the population are adjusted so as to sustain the number of pulsars 
observed at longer periods. 3, Initial periods can be inferred for pulsars that reside within 
composite supernova remnants which are powered by the pulsar spindown energy (van der 
Swaluw & Wu 2001). The derived P via these 3 ways are 10-60 ms, ~ 300 ms, and 37-82 
ms (sometimes several hundred milliseconds), respectively. In addition, Lai et al. (2001) 



1 Pulsars could be neutron stars or strange stars (e.g., Xu 2001). We discuss in this paper both possibilities 
of proto-pulsars being proto-neutron stars and proto-strange stars. 
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studied the implication of the apparent spin-kick alignment in the Crab and Vela pulsars, 
and found that the initial period should be < 1 ms in both the electromagnetic rocket model 
and the hydrodynamic natal model, or < 10 s in the asymmetric neutrino emission model, 
in order to produce a high kick velocity. 

In this paper, we try an alternative effort to derive the initial period through a relation 
between the initial period and the dynamo-originated magnetic field of pulsars. The mag- 
netic field B is assumed not to change in our discussion, since both observation and theory 
imply that a pulsar's 5-field does not decay significantly during the rotation-powered phase 
(Bhattacharya et al. 1992 and, e.g., Xu & Busse 2001). Therefore the field strength observed 
for a normal pulsar represents its initial value. On another hand, the strong magnetic field is 
supposed to be created by dynamo action during the proto-pulsar stage (Thompson & Dun- 
can 1993, Xu & Busse 2001). We may obtain the initial period if we can find an estimated 
relation between B and Pq during the dynamo episode. 

2. The model for calculation of pulsar initial periods 

2.1. The model 

Dynamo action may occur during the first a few seconds of both neutron stars and 
strange stars (Thompson & Duncan 1993, Xu & Busse 2001). In the magnetohydrodynamic 
(MHD) mean field dynamo theory, the helical turbulence (leading to the a effect), differ- 
ential rotation (the Q effect), and the turbulent diffusion (the j3 effect) are combined (e.g., 
Moffatt 1978, Blackman 2002). The dynamo-created magnetic fields are complex and not 
very certain for specific astrophysical processes due to the nonlinearity of the MHD equa- 
tions. Nevertheless it is generally believed (Blackman 2002 and references therein) that a 
seed field begins to grow kinematically until a saturation (energy equipartition between field 
and fluid) phase when the growth slows down due to the back action of field on fluid in a 
small scale (at or below a typical length of turbulent eddies). 

The dynamo-originated field energy in this small scale is provided by both of convection 
and differential rotation. The differential rotation energy is E dii ~ 10 46 P ~ 2 erg (from eqs.(8) 
and (37) of Xu & Busse 2001). For proto-strange stars, the convection energy is E con ~ 
4irR 2 Lpv 2 ~ 10 48 f| erg (Xu & Busse 2001), where R ~ 10 6 cm is the stellar radius, L ~ 10 5 
cm the convective thickness, p ~ 4 x 10 14 g/cm 3 the density of strange quark matter, and the 
convective velocity v — v$ x 10 8 cm/s could be either ~ 10 8 cm/s for a large-scale convection 
or ~ 10 5 cm/s for a local turbulent eddy. For proto-neutron stars, a similar estimate gives 
-Econ ~ 10 48 erg (Thompson & Duncan 1993). We assume that the differential rotation 
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energy dominates, and thus neglects the convective energy, during the episode when dynamo 
action works, since E dif > E con as long as P < 100 ms. The saturation magnetic field in 
small scale, -B sat , satisfies 

Bit _ (1) 



8tt R 3 
Since E di f oc Pq 2 , we come to 

B sat = AP~\ (2) 

where A is a constant for a pulsar-like compact star with typical mass and radius. 

The growth of a much large scale magnetic field (e.g., the observed dipole one, B d ) by 
the a effect only is almost impossible because the magnetic helicity conservation eventually 
leads to a so called "catastrophic" quenching. Nonetheless, differential rotation may play a 
great role to generate large scale fields since it can regenerate a larger scale toroidal field. 
Assuming that the dipole magnetic field, B d , is a constant fraction of the saturation one, we 
have 

B d = aP -\ (3) 

However, it is possible that the ratio of B d to B S3it is a function of P . The Rossby number, 
R Q = Po/r con (the convective overturn time r con ~ 1 ms for both proto-neutron stars and 
proto-strange stars), is an effective dimensionless number to describe the importance of 
rotation, and may also be representative of differential rotation. Therefore the ratio may 
depend on a function of R Q since differential rotation is essential for large sale field generation. 
Observationally, the magnetic activity, which reflects the field strength in the atmospheres, 
of G-M dwarf stars correlate with ^-°- 55 ± a10 (Simon 1990). 2 Considering eq.(2), we can thus 
assume 

B d = a'Po" 1 - 5 . (4) 

Eqs.(3) and (4) are two estimated relations between field and initial period, which will be 
used respectively to derive pulsar initial period in section 2.2. 

After the proto-pulsar phase when dynamo action works, a pulsar may spindown because 
of the braking torques due to magnetodipole radiation and the unipolar generator (Xu & 
Qiao 2001), with a rate 

" = " W tf (5) 

where / ~ 10 45 g cm 2 the momentum of inertia, rj E (1,3) being a function of B d , P, and 
inclination angle a. For the simplicity, we assume i] = 1 in this paper, which is also the case 



2 The energy source for field creation in these dwarf stars could be mainly the convective energy, which 
(and thus B sat ) may be less relevant to the spin period. 
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of conventional magnetic dipole model for pulsars with a = 90°. We then have the magnetic 
field B 12 = B d x 10 12 G by 



The magnetic field is assumed not to change; the field in eq.(6) is thus the same as in eqs.(3) 
and (4), by which the initial period can be obtained. 

We can also have the rotation period P(t) as a function of time t, by integrating eq.(5) 

P(tf = P(0) 2 + (B 2 12 t, (7) 

where P(0) = P , and 

le^io 24 „ 

C = 3^ ^ 1-95 x 10- 15 s. 

Time derivative of eq.(7) leads to (Bf 2 = 2P(t)P(t), we then have the age T, 

p2 _ p2 

T=- t±. (8) 

2PP 

Eq.(7) may be useful to study pulsar "current" in the P — P diagram. 



2.2. Calculations based on pulsar data 

There are three approaches to obtain pulsar ages. The first is the characteristic age 
t = P/2P, which could be derived from eq.(8) if Po <C P. Second, for a pulsar associated 
with a supernova remnant (SNR), the kinematic age of the SNR could be regarded as an 
estimation to the pulsar age. Third, when an SNR is identified to be the remnant of a 
historically recorded supernova, the age of the pulsar within the SNR is then precisely 
known. 

Among the youngest pulsars which locate in SNRs (see table 2 of Camilo et al. 2002), 
only the Crab pulsar is confidently believed to be associated with SN1054, the other two, PSR 
J1181-1926 (in the SNR Gll.2-0.3) and PSR J0205+6449 (in the SNR 3C58) are suggested 
to be probably associated with SN1181 and possibly with SN386, respectively (Clark & 
Stephenson 1976). Therefore, only the real age of the Crab pulsar is exactly known in all 
the discovered pulsars. 

For the Crab pulsar, with P = 33.4 ms, P = 4.2 x 10" 13 s/s and T = 948 yr, its initial 
period is then P = 16.8 ms through eq.(8). Calibrated by the Crab pulsar, eq.(3) can then 
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be rewritten as 

B 12 = 0.0637P -\ (9) 

which provides a way to calculate the initial period. The ATNF pulsar catalogue has released 
~1,300 radio pulsars with P and P measured (http:/ /www. atnf.csiro.au/research/pulsar/catalogue). 
In our calculation, millisecond pulsars are excluded, for their magnetic fields may decay sig- 
nificantly during the recycling phase. In this paper, the criteria used to exclude a millisecond 
pulsar is (a) P < 60 ms, (b) T < 10 9 years, and (c) B ~ 10 8 ~ 9 G. Besides millisecond pul- 
sars, those pulsars in binary systems and those having negative values of P are also excluded. 
Totally 1,069 pulsars are selected for our calculation. 

The initial periods are derived via eqs.(6) and (9). The result is presented by the 
histograms for P (left panel of Fig.l). It is found that 60% pulsars have initial period less 
than 60 ms, exhibiting an peak at the interval 20-30 ms, as shown in Fig.l. It is worth 
noting that data of 997 pulsars are shown in the histogram. For the other 72 pulsars, the 
calculated Po is greater than P, hence results in a negative value of T, which means that 
eq.(9) may not be applicable to these pulsars. 

Similarly we can also find the coefficient a' = 8.24 x 10 -3 in eq.(4) by the Crab pulsar 
calibration. Employing this relation, the initial periods are calculated, which is presented 
in the right panels of Fig.l. There are 23 pulsars which have negative calculated T, so that 
totally 1046 pulsars are plotted in the histograms. The distribution peak of P is also between 
20 and 30 ms, but more pulsars tend to have smaller P . Statistically, the distributions of 
characteristic age r and the age T do not show significant difference. This is not surprising 
since most pulsars have their initial period Pq <C P. 



3. Conclusion and discussion 

The statistical distributions of pulsar initial periods are inferred based on the estimated 
relation between the initial periods and the magnetic dipole magnetic fields. It is found that 
proto-pulsars are very likely to have rotation periods between 20 and 30 ms, and that most 
of the pulsars rotate initially at a period < 60 ms. The pulsar initial periods derived in this 
way are comparable with the previous estimates by the first and the third methods, but are 
significantly smaller than that by the second method (section 1). 

When looking at the pulsar distribution in the P — P diagram, one generally concludes 
that the magnetic fields decays at a time scale of 10 6-7 years since the number density of 
pulsars with low fields and old ages are larger than that of pulsars with high fields and 
young ages. However, calculations of ohmic decay of dipolar magnetic fields by Sang & 
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Chanmugam (1987) suggest that neutron star magnetic fields may not decay significantly 
during the rotation-powered phases. The field of a strange star can also hardly decay because 
of the high magnetic Reynolds number and possible color superconductivity (e.g., Xu & Busse 
2001). Observations indicate that the field decays only if the pulsar is in the accretion phase 
(e.g., Taam & van den Huevel 1986). If the field does not change significantly, how can 
we understand the apparent "field decay" in the P — P diagram? We may explain this 
discrepancy by studying pulsar current in the P — P diagram, according to eq.(7) with the 
inclusion of rj ^ 1 in different emission models. Certainly the braking index varies when a 
pulsar evolves in this way. Further studies on this topic will be necessary and interesting. 

There are a few pulsars which have initial periods < 5 ms in our calculation; for instance, 
5 pulsars with P < 5 ms appear in the right histogram of Fig.l based on eq.(4). This result 
can not be understood in the conventional opinion of r— mode instability, 3 if it is not caused 
by statistical error or other factors being inherent in our P Q calculations. However, as 
addressed in Xu & Busse (2001), the r— mode instability may be inhibited by differential 
rotation and strong magnetic field; the problem could thus be solved. 
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Fig. 1. — Histograms for the initial period P . Left: P is calculated from B u = 0.0637P _1 , 
in which data of 997 pulsars are used and each bin represents an interval of 10 ms. Right: 
Pq is calculated from Bu = 8.24 x 10~ L5 P -1 ' 5 , in which data of 1046 pulsars are plotted and 
each bin represents 5 ms. 



